Nonequilibrium dynamics of polymer translocation and straightening 
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When a flexible polymer is sucked into a localized small hole, the chain can initially respond 
only locally and the sequential nonequilibrium processes follow in line with the propagation of the 
tensile force along the chain backbone. We analyze this dynamical process by taking the nonuniform 
stretching of the polymer into account both with and without hydrodynamics interactions. Earlier 
conjectures on the absorption time are criticized and new formulae are proposed together with time 
evolutions of relevant dynamical variables. 
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I. INTRODUCTION 



A long flexible polymer is one of the representative 
examples of soft matter. A common feature of soft mat- 
ter is a presence of mesoscopic length scales, which is, 
in many respects, responsible for their unique proper- 
ties such as the high susceptibility. For dilute poly- 
mer solutions (with each chain made of N succession of 
monomers of size b), this corresponds to the Flory radius 
i?o = bNy of individual coils, which serves as a basis 
for the scaling theory From the estimation of the 
elastic modulus ~ ksT / one can realize an important 
consequence that a long chain is readily exposed to sig- 
nificant distortions, such as stretching and compression, 
by rather weak perturbations. 

Although the extension of a polymer in various flow 
fields or by mechanical stretching have been extensively 
studied [J 0, [H, 0, H, @, most attentions so far have 
been focused on equilibrium or steady state properties. 
One can also ask the dynamical process from one steady 
state to the other induced by sudden change of external 
field Q, which would be important in relation to recent 
development of micromanipulation techniques. For in- 
stance, imagine that an initially relaxed polymer is sud- 
denly started to be pulled by its one end (see Appen- 
dex [Cj) . If the force is sufficiently weak fn < 1, a chain 
as a whole follows at the average velocity v ~ f/(rj s Ro) 
(j] s is the solvent viscosity) with keeping the equilibrium 
conformation. (Here and in what follows, we denote the 
dimensionless force f x = fx/iksT), where x has the di- 
mension of the length.) For large force fn > 1, how- 
ever, only a part of the chain can respond immediately, 
while remaining rear part does not feel the force yet. As 
time goes on, the tension propagates along the chain, 
which alters chain conformation progressively, and the 
steady state is reached after a characteristic time. At 
room temperature, the critical force / ~ (fcsT)/ Rq is es- 
timated to be on the order of pico Newton for a flexible 
chain with N = 100 ~ 1000, comparable to the usual 



force range in the single molecule manipulation experi- 
ments with atomic force microscopy and optical tweezer. 
The force generated by molecular motors also falls into 
this range. This implies a possible importance of such 
a nonequilibrium response in many biological as well as 
technological situations. 

In the present paper, we illustrate such a nonequilib- 
rium response using an example of polymer absorption or 
aspiration into a small spot. Our target here is the dy- 
namical process, in which a polymer is sucked into a lo- 
calized hole. This is different from the pulling the chain's 
one end Q , but similar in the way how the chain responds 
to the local force, and, in fact, relevant to the dynamics 
of polymer translocation through hole [§, [l(| and the ad- 
sorption process to a small particle. (In this case, the 
force / is related to the chemical potential change A/i 
due to the absorption via / = A/x/6.) Although the phe- 
nomenon of polymer translocation has been an active re- 
search topic in the past decade as a model for biopolymer 
transport through a pore in membrane, our current un- 
derstanding for the dynamics, in particular the strongly 
driven case, is restrictive. So far, the scaling estimates 
of the characteristic time r for absorption process in im- 
mobile solvents have been proposed (g, [ll[. Grosberg et. 
al. argued the absorption time for a Rouse chain on the 
ground that Rouse time tr = t Nq is the solitary rele- 
vant time scale (to — i] s b 3 / (ksT) is a microscopic time 
scale), and all other relevant parameters appear in the 
dimensionless combination /r , thus, 



TR(j>(fR ) 



(1) 



The scaling function is determined from the require- 
ment that the speed of the process, Nq /r must be linear 
in the applied force, leading to 



Vsb 3/2 

~T 



(2) 



This result was interpreted as a sequential straightening 
of "folds" [III- For a chain with excluded volume studied 
by Kantor and Karder, a similar scaling argument leads 
to 1 
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(3) 
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where Flory exponent v = 1/2 for a chain in 9 solvent, 
while in good solvent v ~ 3/5 (in space dimension T> = 3) 
and = 3/4 (in V = 2). 

We tackle this problem through different approach by 
explicitly considering the dynamics of the tension propa- 
gation. This enables us to unveils the physics behind the 
nonequilibrium driven absorption and go beyond the pre- 
vious works by taking the excluded volume effect and/or 
hydrodynamics interactions into account. We indeed find 
that eq. (I2J) , d3j) are not generally correct as a conse- 
quence of the salient feature of flexible chain that the re- 
sponse to the aspiration or stretching force is nonuniform 
both in space and time. At very strong driving, the chain 
finite extensibility matters, and this leads us to propose 
three distinguished regimes for the absorption process de- 
pending on the degree of forcing. Interestingly, we shall 
see that eq. ©, ([3]) are recovered in the limit of very 
strong forcing only. In addition to the absorption time, 
we can also predict the time evolutions of dynamical vari- 
ables governing this nonequilibrium process. Below, the 
problem of the dynamical response is formulated with 
basic equations and the absorption dynamics is analyzed 
in Sec HTAl Then, Sec. liTBl and liTCl are devoted to dis- 
cussions on the effect of the finite chain extensibility and 
hydrodynamic interactions, respectively. Summary and 
future perspectives are given in Sec. IIIII Some technical 
details and mathematics are given in Appendix lAl and [Bl 
Another example of nonequilibrium response, i.e., a sud- 
den pulling of a chain by its one end, is briefly discussed 
in Appendix [Cl 



II. FORMULATION 

A. Dynamical response to strong forcing 

Now, let us imagine the moment, when one end of the 
chain arrives at the attractive hole at the origin (Fig. [1} . 
The fact that the first monomer is strongly pulled in- 
fluences the rear vicinity immediately, but not far away. 
If we notice the subunit consisted of first <?o monomers 
of size yo = bg^, this subunit starts to move with the 
average velocity vq ~ f /(r] s yo), provided that the defor- 
mation is insignificant in this scale. The size of such a 
subunit is deduced by comparing the longest relaxation 
time of the subunit Togo" to the velocity gradient vq / yo 



9o 



k B T 
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(4) 



and this constitutes the initial condition. 

The absorption process proceeds with time; at time 
t, the tension is transmitted up to iV(t)-th monomer, 
while M(t) front monomers are already absorbed (Fig. 
1). This indicates that one can regard that the chain 
portion at (—R(t) < x < 0) takes a steady state confor- 
mation moving with the average velocity v(t). Here the 
drag force builds up along the chain starting from a free 




FIG. 1: A dynamical response of a polymer initially at rest 
to the strong aspiration. A chain is pulled from (x < 0) re- 
gion to the aspiration spot located at the origin. Monomers 
in gray region are under the influence of the tensile force, 
while monomers in the rear part is yet relaxed. The distance 
between the boundary (associated with tension propagation) 
and the origin and the size of the largest blob at the bound- 
ary are denoted as R(t) and yn(i), respectively. Monomers 
are sequentially labelled as 1, 2, ■ ■ ■ No from one end (arriving 
at the origin initially) to the other. M(t) and N(t) are, re- 
spectively, labels of the monomers residing at the origin and 
the boundary at time t. (top) Nq" < fb < 1; (bottom) 
1 < < Nq, in which the anterior part (— l(t) < x < 0) is 
completely stretched. 



boundary (x = —R(t)) to the origin, which makes the 
overall chain shape reminiscent to a "trumpet" ||. It 
means that the length scale y(x) set by the large tensile 
force, above which the chain is substantially elongated, 
is position dependent, thus, the elastic behaviour of the 
chain is described as a sequence of blobs with growing 
size y(x) — bg(x) v . This leads to the following local 
force balance equation; 



y{x) 



r) s v(t) \x + R(t) 



-R(t)+y R (t)<x<0](5) 



where we introduce the cut-off length yn(t), which signi- 
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fies the size of the largest blob at the free boundary; 



y{-R{t)+y R {t))~y R {t) 



y R (t) 



k B T 
rj s v(t) 



1/2 



[at x = -R(t) + y R (tJ$) 



(Note that throughout the paper, we neglect the loga- 
rithmic factor associated with the friction of asymmetric 
objects in low Reynolds number [l2j regime as well as 
other numerical coefficients of order unity unless speci- 
fied.) The mass conservation reads 



R(t) 



p(x)S(x)dx + M(t) = N(t) 



(7) 



where p(x) ~ g(x)/y(x) 3 is the monomer density (y(x) = 
bg{x) v ) and S(x) ~ y{x) 2 is the cross-sectional area of 
the conformation. Equation ([5]) and ([7]) constitutes ba- 
sic equations supplemented with the following statistical 
relation available from the conformation at rest (t < 0) 



bN{tf = R(t) 



(8) 



and "boundary conditions" both at the free boundary eq. 
([6]) and the origin 



r) s R(t)v(t) ~ / 



(9) 



These conditions express the force balance at the free end 
and the total force balance, respectively. 

After casting Eq. ([7|) in the differential form 



[pS\x= 



(dR{t) 
R(t)) — 7T- 



V dt 



v{t) 



dN(R) 
dt ' 



(10) 



we obtain the following equation for the tension propa- 
gation (see Appendix); 



(R(t)\^ $(R(t)) = t 
\ Ro ) $(i? ) n 



(11) 



With the symbol x possessing the dimension of length 
(recall once again the definition of the dimensionless force 
f x ), the function <j>(x) is defined as 

S(a:) = l-B Q (f x )^ + (B -l)(/*r^ i (12) 
* l-B (f x )^ (13) 

where -Bo is the numerical coefficient of order unity and n 
is the time, when the other side of the chain end reaches 
the steady state and gets set into motion due to the ten- 
sile force; i?(n) = Rq, thus, 



n = T Z (fRo) 2 » v ®(Ro) 



(14) 



where Tz — t Nq 1 ' is the Zimm time. We notice that the 
time Ti obtained here satisfies the scaling form of eq. (H|) 
with the replacement of tr by tz- 
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FIG. 2: Time evolutions of number of absorbed monomers 
M(t) at various forces for the chain length iVo = 10 3 under 
good solvent condition (v = 3/5). Inflection-like points mark 
the end of the tension propagation stage (t = n and M = 
Mi). 



The number of the monomers absorbed is calculated 



as 



M(t) = / [pS] x = v(s)ds 
Jo 



(15) 



which can be shown to reach Mi = N {1 - {}r ) {v ~ 1),u ) 
at t = T\. After that, the whole yet non-absorbed part 
of chain is under the tension and pulled toward the hole, 
thus, the evolution of M(t) is governed by 



(16) 



where L(t) is the long axis length of the non- absorbed 
chain. This leads to UM 



M(t) = M 1 + f Ra (fb)-- 



1 " 1 " fRo 



t-n 

TZ 



[17) 



[t > n ] 



The absorption process completes at time r = t± + T2 
with 



r 2 ~ T Z (fR ) 



(18) 



The time evolution of M(t) is plotted in Fig. [2] for 
iVo = 10 3 under various forces. One can see that the 
first stage of the tension propagation dominates the most 
process under the large force f Ro > 1. In this case, from 
eq. (|14j) . the absorption time is approximated as 



T — Tz UrS 2 ^ 



TO (fb 



-N~ 



(19) 



In more specific form, r ~ TV 7 / 5 / J 2 / 3 for a chain in good 
solvent [u = 3/5) and r - iV 5 / 4 // 1 /2 for a chain in q 
solvent, which are apparently different from earlier con- 
jectures (eq. ©, (O). As eq. (TTTj) indicates, time evolu- 
tions of dynamical variables with suitable normalization 



4 




FIG. 3: Normalized time evolutions (double logarithmic scale) 
of number of absorbed monomers M(t/n)/Mi at various 
scaled forces in the tension propagation stage under good sol- 
vent condition (y — 3/5). 



are represented by master curves parametrized by values 
of fn in this tension propagation stage (thus, almost all 
the absorption process). Figure [3] exemplifies the time 
evolution of M(t) with the normalization M(t)/M\ and 
t/n, where M x = N (l - (f Ra ) {u ~ 1),u ). The larger the 
scaled force fn , the sooner the master curve approaches 
the asymptotic line with the slope (1 + i/)/(l + 3^). 



B. Effect of the finite chain extensibility 

The above argument is valid as long as /& < 1. For 
larger forces, as one can see from eq. (HJ), the chain close 
to the origin is completely stretched (see the bottom of 
Fig. [I}. To include such an effect, let us begin with the 
limit of strong force (f b > Ng)\xft. Then, at any mo- 
ment, the entropic coiling is completely irrelevant for the 
tensed part of the chain. We can repeat the same ar- 
gument by setting y(x) — b for (— R(t) < x < 0), and 
find 



with 



n = r (fb) 



T R (fRoY 



(20) 



(21) 



where Rouse time tr = T N (R /b) 2 (for a chain with 
Flory exponent v) appears as a characteristic time re- 
flecting the fact that the backflow effect is nearly negligi- 
ble in this strong pulling limit. It is noted that eq. (12"T]) 
coincides with eq. ^ and ([3]). 

For the intermediate case of 1 < fb < Nq , this process 
of strong pulling limit is valid until the time r*. This 
time is determined by monitoring the monomer at the 
free boundary; the velocity of the tense string of length 



l(t) = N(t)b is given by 
dN(R) 



v(t) 



dt 




(22) 



where we have utilized eq. ([5]) and (|2"D|) . This quantity 
decreases with time and at t = r* , the drag force for the 
free boundary monomer becomes comparable to ksT /b. 



Vsbv(r* 



From this, we find 



= r (A)- 



(23) 



(24) 



After r*, the conformation of the chain under tension 
is a completely stretched string of length l(t) followed 
by coiled subunits with growing size [6(. For the latter 
part, we can apply the aforementioned analysis for the 
trumpet, where the "boundary condition" is imposed not 
at the origin but at x — l(t). Setting R(t) = l(t) and 
y(t) — b in the local force balance equation (eq. (HJ)) 
yields 



l(t) ~ R(t) 



Vsbv(t) 



R(t)(i - (hr 1 ) 



(25) 



where the second equality utilizes the total force balance 
r] s v(t)R(t) ~ /. More precisely, eq. (|25|) should be writ- 
ten as 



P 



Kt) 



R(t) 



1 - c (fbY 



(26) 



where Cq is the logarithmic factor dependent on l(t) and 
R(t). Ignoring this small correction allows us to derive 
the following form of the tension propagation dynamics 
at t > r* 



Tl 



&*(R(t)) 
$*(Ro) 



(27) 



where the function <fr*(x) and the time r* are respectively 
defined as 



$*(x) 



T l - Tl 



l—p 2- 

T Z (fRoY 



— B (f x 
■**(i2o) 



1-P Z 



(28) 
(29) 



Absorption time is approximated by the time t\ when 
the tension reaches the other end, thus, from eq. 
and 



tzUroY 



(30) 



Utilizing the asymptotic behaviours of eq. (|26p . i.e., p — > 
for fb — » 1, and p — * 1 for fb >> 1, eq. (|3^|) shows 
smooth crossovers to the trumpet regime eq. (fl"9| at 



fb — > 1 and strong pulling limit eq. (f2"Tj) at fb 



NX. 
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C. Effect of hydrodynamic interactions 

At certain situations, hydrodynamic interactions may 
become screened, i.e., a polymer confined in narrow slit 
or in melt of short chains. Then, a question arises; what 
is the effect of the hydrodynamics in this dynamical pro- 
cess? To answer this, let us "switch off' the induced 
flow of solvents. Then, the only requisite alteration is 
the dissipation mechanism, which becomes local and in- 
dependent of the chain conformation. The local force 
balance is, instead of eq. ([5|), written as 



rj s bv(t) 



where g(x) is related to y(x) through y(x) = bg(x) u . 
largest blob size at the free boundary is 



VR(t) 



k B T 
b 2 rj s v(t) 



(31) 



The 



(32) 



One can analyze along the same line as the chain in mo- 
bile solvents. In partucular, t± for the moderate forcing 
(fb < 1) is obtained as 



n 



(33) 



with 



(35) 



which coincides with the scaling form of eq. |T]) and 
should be contrasted with eq. (|14| with hydrodynamic 
interactions. On the other hand, the result in the limit 
of the strong forcing (/^ > Nq ) is not altered, and t\ is 
given by eq. (|2ip ■ It is in this limit only that the require- 
ment t ~ / _1 is fulfilled reflecting the saturation of chain 
deformation due to the complete stretching, therefore, 
the earlier conjectures are approved [1, [ll[ . For smaller 
forces, the soft elasticity of the chain results in more in- 
volved responses as we have seen; the stronger driving 
results in the more intense chain deformation, and this 
deformation behaviour affects the absorption dynamics. 



III. SUMMARY AND PERSPECTIVES 

There would be many practical situations, in which 
the externally imposed velocity gradient exceeds the in- 
verse relaxation time of long polymers. If the external 
field acts locally, effects associated with nonequilibrium 
response, i.e., the propagation of the tensile force along 
the chain, are expected to show up. We focused on the 
problem of absorption or aspiration into a localized hole 
and demonstrated how such a process can be physically 
described. 

One of the most relevant situations of aspiration dy- 
namics studied here is the polymer translocation through 



a pore. We should note, however, in the problem of 
biopolymer transport through a membrane pore, the role 
of specific interactions may become essentially impor- 
tant [§, [l^, EH • We neglect all the complications associ- 
ated with such a factor and focused on universal aspects 
as a consequence of a polymeric nature, in this sense, the 
present analysis may be regarded as an "ideal" version in 
view of the relation with the problem of polymer translo- 
cation. Such an "ideal" situation would be now experi- 
mentally feasible thanks to the advance in nanoscale fab- 
rications fLU ]. 

The advantage of the present framework is rather wide 
range of applicability to related problems. As mentioned 
in Sec. [H if the chain is suddenly pulled its one end, the 
response is nonuniform both in space and time. The re- 
sultant transient dynamics of the chain extension can be 
analyzed in a similar way (see Appendex |C|) . The same 
physics is also expected in the escape process of a con- 
fined polymer from a planner slit [15j . Considering hier- 
archical structures common in polymeric systems, such 
nonequilibrium dynamics in a single chain level would be 
expected to show up in macroscopic material properties, 
too. We hope that the present study provides the basic 
insight involved in the driven nonequilibrium process of 
polymer absorption and its related problems and future 
inve stig ations including the comparison with the numer- 
ical [18[ and even real experiments would be valuable. 
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APPENDIX A: CONTINUITY EQUATION 

In this appendix, we shall discuss the relation between 
integral and differential forms of the mass conservation 
equation. We define h(x) — p(x)S(x) and J h{x)dx = 
H{x) for the concise notation. Then, the integral form 
of the mass conservation (eq. ([7])) is written as 



H{0) - H(-R(t)) + M(t) = N(t) 



(Al) 



The variation of this equation with the variable transfor- 
mation from x to u = x — v(t)t leads to 



dH \u=-v(t)t - 

h(0) x [-vdt] 
dN 



dH 



u=-R {t) ~ vi t )t + dM (t) = dN(t) (A2) 
h(-R(t)) x [d(-R(t)) - v(t)dt] + dM(t) 



dR 



dR(t) 



(A3) 



Since h(0) x v(t)dt is the number of absorbed monomer 
during the time interval dt, the first and the last terms in 
left hand side are canceled out. The resultant equation is 
the differential form of the mass conservation (eq. (JTUJ) ) . 
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APPENDIX B: SOLVING DIFFERENTIAL 
EQUATIONS 

In this appendix, we illustrate some technical de- 
tails how to solve the differential equations. Let us see 
eq. ([TO]), which is coupled with eq. @, © and j9|). Using 
the definition of density p(x) and the cross-sectional area 
S(x) (given below eq. 0) and also eq. ([8]), one obtains 



(i?(t)^ -{VRit)) 1 ^) 



dR(t) 
dt 



By putting force balance conditions (eq 
eq. (|Bip is transformed to 



(y R ( t) )^v(t) (Bi) 

I and ©), 



R(t)^~ - 



k B T 
f 

k B T 



f 



R(t) 
R(t) 



dR(t) 
dt 

f 

T) s R(t) 



(B2) 



The solution of eq. (|B2[) with the initial condition eq. ^ 
is 



R(t) 



R(t) 



J Vo 



-Both)- 1 



2] «(*) 



J Vo 



t 



(B3) 



where we introduce the numerical coefficient Bo of order 
unity to replace the relation symbol ~ with =. Using the 
function <f> (eq. (|12[) ). the above equation is rewrriten in 
the compact form 



(B4) 



The time t\ (eq. (fl4"l) ). when the tensile force reaches 
the chain end is obtained by putting R(t) = Rq in this 
equation. Equation is obtained after normalizing 
t by T\. Equations (f2"T)) and (|30[) in the case of strong 
pulling (1 < fb < Nq) are obtained in the same way just 
by replacing the lower limit of the integral in eq. (|B3[) 
with l(t) (eq. (|2l))0 and the right hand side with (t — 



APPENDIX C: PULLING ONE END 

Here, we shall briefly discuss another example of 
nonequilibrium response, i.e., the transient dynamics 
of the chain stretching pulled by its one end (Fig. d]). 
The difference between the absorption and the present 
stretching process was pointed out by Kantor and Kar- 
dar [f| , which is clearly recognized by comparing Fig. [JJ 
and Fig. in the former, the site of action is fixed in 




FIG. 4: A dynamical response of a polymer initially at rest 
to the strong pulling by its one end (Nq 1 ' < f b < 1). A coor- 
dinate is defined in such a way that the origin is the position 
of one chain end at time t = and that end is started to be 
pulled to right (x > 0) direction. Monomers in gray region 
are under the influence of the tensile force, while monomers in 
the rear part is yet relaxed. The distance between the bound- 
ary (front of the tension propagation) and the origin and the 
size of the largest blob at the boundary are denoted as R(t) 
and yn(t), respectively. Monomers are sequentially labelled 
as 1, 2, • • • No from the pulled end to the other. N(t) is the 
label of the monomer residing at the boundary at time t. 



space and the chain portion after crossing this point gets 
relaxed, while in the latter case, the site of action is 
moving with the velocity v(t) = dl(t)/dt and the chain 
portion after crossing the origin is still tensed and con- 
tributes the friction. 

In this case, too, one can distinguish three different 
regimes depending on the applied force; (i) trumpet 
(fb < 1 < Iro), (h) intermediate (1 < f b < N ) and 
(iii) strong limit with complete stretching (/& > No) by 
monitoring the force acting on the last monomer at the 
free boundary (Note that the border between regime (ii) 
and (iii) is different from the absorption case). Just as 
the case in the absorption dynamics, one can write down 
basic equations. The local force balance; 



y(x) 



k B T ( 1 



j] s v(t) \x + R(t) 
The mass conservation; 

W 



-R{t)+y R (t)<x<l(1$pi) 



R(t) 



p(x)S(x)dx = N(t) 



(C2) 



One also needs eq. ((8]), the condition for the free bound- 
ary (eq. ((U) and the total force balance; 



y(x = l(t)) 

& l(t)+R(t)~ 



k B T 



T) s v(t) 



(C3) 
(C4) 



By substituting eq. (|C1[) into eq. (|C2[) and using eq. Q 
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and (|C4|) . one obtains the relation between v(t) and N(t); 

2v-l 



f 



k R T 



Tjv(t) 



bN(ty 



(C5) 



The steady state velocity v ss after the arrival of the ten- 
sile force at the other end is found by setting N(t) = No 
in this equation; 



VsR 2 



(C6) 



On the other hand, there is another relation between 
v(t) and N(t), which is available from eq. (JSj) and (|C4jl ; 



bN(tf 



f 



rj s v(t) 



-w 



f - [ v(t')dt' 



T) s v(t) 



(C7) 
(C8) 



Combining cq. (IC5|) and (|C8j) leads to an integral equa- 
tion for v(t), solution of which is 



/ 



r] s v(ty 



1-Bo 



rj s k B Tv(t) 



1-1/ 



(C9) 



(where we dropped a small term associated with the ini- 
tial condition) . The time t\ for attaining the steady state 
is found by setting v(t) = v ss in this equation; 



n - Tz(fRo) - l - B {f Ro ) 



(CIO) 



After rescaling the time and velocity by t\ and v ss re- 
spectively, eq. (|C9[) is rewritten as 



'«(*) 
. v ss 



- (CH) 
n 



Above argument is valid when the applied force does 
not exceed the threshold /& < 1. For larger forces, the 
nonlinear effect associated with the finite chain extensi- 
bility becomes apparent. In the limit of strong forcing 
fb > No, the lateral chain size of the tensed portion be- 
comes just y(x) — b and the analysis becomes very easy 
as in the case of absorption dynamics. Repeating the 
same argument as above, one finds 



_ k B T fb 
VsbN ~ Aq 



,6 2 A n 2 



n 



= T R (f Ro )- 2 f b N 



(C12) 



(C13) 



Vss 
V(t) 



I -Bo 



v ss N 
v(t) 



n 



(C14) 



The asymptotic form of the velocity v(t) of the pulled 
end is 



(C15) 



In the intermediate case 1 < fb < No, one should have 
crossover between above two regimes just like the absorp- 
tion dynamics. Equation (|C13[) was proposed and con- 
firmed by Monte Carlo simulation in the limit of strong 
forcing 8]. It is important to notice that in the case of 
pulling one end, too, the earlier conjecture is approved 
for the strong force limit only. 
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